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Summary 

A study has been conducted to determine how well the USM3D unstructured Euler solver can be 
utilized to predict the flow over the High Speed Research Reference H configuration with an ultimate goal 
of prediction of sting interference so after body closure effects may be evaluated. This study has shown 
that the code can be used to predict the interference effects of a lower mounted blade sting with a high 
degree of confidence. It has been shown that wing and fuselage pressures, both levels and trends, can be 
predicted well. Force and moment levels are not predicted well but experimental trends are predicted. 
Based upon this, predicted force and moment increments are assumed to be predicted accurately. 
Deflection of the horizontal tail was found to cause a non-linear increment from the non-deflected sting 
interference effects. 

Introduction 

The pursuit of interference free wind tunnel data has kept many an engineer gainfully employed since 
the wind tunnel was invented. Corrections have been developed to take into account the effect of wind 
turinel walls and for blockage. The effect of the sting in configuration aerodynamics is minimized with the 
use of a straight, circular sting that intersects the aft end of the fuselage. This arrangement results in little 
interference, especially at lower angles-of-attack, since the sting can be treated as a streamtube. 
However, the aft body of the configuration being tested sometimes requires significant modification to 
accommodate this type of sting leading to the requirement for aft-body corrections. One method of wind 
tunnel model support that allows for the true after body to be modeled is a blade support that can intersect 
either the upper or lower fuselage surface forward of the empennage. This support structure, though, can 
no longer be treated as a streamtube and may cause significant interference on the configuration 
aerodynamics, especially at yaw angles and transonic speeds. 

The purpose of this study is to first, validate the USM3D unstructured Euler solver for the prediction of 
flows about the High Speed Research Reference H configuration. A complete configuration will be 


analyzed, including engine nacelles, horizontal tail and vertical tail. The second objective of this study is to 
numerically determine the affect of a lower, aft-mounted blade sting on the configuration aerodynamics. 
Forthis first study, sting effects will be determined only for the non-yawed configuration. These increments 
will be used to correct experimental data for sting interference so an empennage effect can be determined. 
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Subscripts: 

HT 
LE 
TE 

Method Description 

The computational method utilized for this study was the USM3D Euler solver developed by Frink, 
Parikh and Pirzadeh 1 . This code utilizes an upwind numerical scheme to solve the time-dependent, three 
dimensional Euler equations on an unstructured tetrahedral mesh. A cell-centered finite-volume 
formulation is employed with both Roe’s flux-difference splitting and Van Leer’s flux-vector splitting 
available for spatial differencing. 2 Solutions are advanced in time to steady-state conditions by means of 
an implicit Gauss-Seidel algorithm described in reference 2. A first order interpolation scheme as well as a 
higher order scheme are also available. 


span, inches 
drag coefficient 
lift coefficient 
moment coefficient 
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fuselage station, inches 
freestream Mach number 
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drag coefficient increment 
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Generation of the unstructured meshes used in this study was accomplished via the GridTool 3 and the 
VGRID 4,5 6 software. These codes develope the computational grid using the advancing front 
technique. Details of each code can be found in the respective references. 

Grid Generation 

Semi-span unstructured grids were constructed for this study about the High Speed Research (HSR) 
Reference H configuration for various leading and trailing edge flap configurations as listed in table 1 . 
Grids were developed for the configuration with a lower mounted blade sting, a configuration with a circular 
aft body and circular sting, as well as for the configuration with no wind tunnel support structure modeled 
(referred to as the “free air" configuration). In all cases the horizontal and vertical tail as well as the engine 
nacelles and diverters were modeled. Computational meshes were also developed about the 
configuration with a -5' horizontal tail deflection. No grids were constructed for a rudder deflection. A 
view of a typical surface mesh is shown in figure 1 . A detailed view of the wing lower surface showing the 
30* leading-edge and 20* trailing-edge flap deflection as well as the nacelles and diverters is presented in 
figure 2. 

Results 

Code Validation 

Solutions have been obtained for the HSR Reference H configuration at the conditions listed in table 1 . 
The locations of experimental wing and fuselage surface pressure coefficients are shown in figures 3 and 
4. Pressure coefficients obtained from the USM3D Euler calculations are compared with unpublished 
experimental data in figures 5 - 31 . (The experimental data was obtained in the National Transonic Facility 
(NTF) at NASA Langley Research Center.) All solutions were obtained utilizing the higher order 
interpolation scheme and flux difference splitting with the exception of the M=0.95, a=8* solutions. These 
solutions (with and without the blade sting modeled) were only obtainable using first order interpolation 
and flux vector splitting. 

The wing surface pressure coefficient plots exhibited in figures 5-9 show comparison of experimental 
data and computational results for an angle-of-attack variation at 0.95 Mach number. The theory predicts 
the experimental data well for all but the higher angle-of-attack cases. Here, as expected, the inviscid 
theory does not predict the apparent leading-edge vortex, a viscous phenomenon, evident in the 
experimental data (figures 8 and 9, Fuselage Stations 46.5 and 49.55). The code shows equally as good 
agreement with experiment for a deflected horizontal tail geometry (figures 10 and 11). Undeflected 
horizontal tail experimental data is used for comparison since deflection of the horizontal tail should not 
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significantly affect the wing surface pressures. At subsonic conditions with large leading and trailing edge 
flap deflections the theoretical method predicts the experimental results reasonably well (figure 12). 

The plots presented in figures 13-17 show good comparisons of predicted fuselage surface pressure 
coefficients with experimental levels for a circumferential distribution of pressure taps at Fuselage Station 
65.306 as angle-of-attack is increased for 0.95 Mach number. Note that the computational method is 
predicting the break in the experimental data in figures 13 - 17 between azimuth angles of 105* and 150*. 
Good comparisons are also obtained for a horizontal tail deflection at transonic conditions (figures 1 8 and 
19). (The experimental data used in these comparisons was obtained for a geometry with deflected 
horizontal tail.) Figure 20 shows that a large trailing-edge deflection at subsonic speeds does not lower 
the quality of the comparisons. In fact, the comparisons at this condition is even more favorable than the 
transonic comparisons with undeflected inboard flaps. 

Surface pressure comparisons along a line 45' below the horizontal extending from F.S. 65.306 to F.S. 
81.04 for angle-of-attack varying from 0’ to 8’ is shown in figures 21 to 25 at 0.95 Mach number. In 
general, the theory agrees well with experiment, predicting the trend of the data well. The effect of a 
deflected horizontal tail is also predicted (figures 26 and 27). There is excellent agreement of computed 
results with experimental data for the subsonic case with large flap deflections (figure 28). These plots 
show there is a significant effect of blade sting interference on the aft-fuselage surface pressures and that 
when the geometry analyzed is consistent with the wind tunnel model, much better agreement with 
experimental levels are obtained. 

Comparisons of predicted lift, drag and pitching moment with unpublished experimental values are 
shown in figures 29 - 30. The experimental data has not been corrected for aeroelastic effects. The Euler 
method predicts higher lift coefficient levels and more negative pitching moments than observed in the 
experiment (figure 29). These differences are consistent with those resulting from small aeroelastic 
effects. However, the trends of both lift and pitching moment with angle-of-attack are captured in the 
computed results. An increase in drag is also predicted at each angle-of-attack over the experimental level 
as would be expected with an over prediction of lift. (Note that the equivalent flat plate skin friction drag 
coefficient for eighty million Reynolds number has been added to the predicted drag coefficients.) 

However, when drag coefficient is plotted against lift coefficient (figure 30), the comparison of predicted 
values with experiment improves significantly. 

These comparisons presented in figures 5 - 30 shows that the USM3D Euler method, while not 
predicting absolute levels, predicts overall trends in the experimental data well. This suggests that this 
code can be utilized to predict force and moment increments due to the blade sting support structure. 
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Sting Interference Predictions 

Fuselage incremental surface pressure coefficients resulting from the interference of the blade sting in 
the computational analysis at 0.95 Mach number are presented in figure 31 (streamwise) and figure 32 
(circumferential) with force and moment increments shown in figure 33. For both the calculations with and 
without the sting, the flaps on only the outboard portion of the wing were deflected -- 10’ down along the 
leading-edge and 3' down along the trailing-edge. The predicted increments are calculated by subtracting 
the levels obtained from the analysis of the configuration without any sting support structure modeled from 
the levels resulting from analysis of the configuration with the blade sting modeled. 

Pressure increments presented in figures 31 and 32 reveal that the effect of the blade sting support 
structure on the aft fuselage lower surface pressures is moderately independent of angle-of-attack. The 
effect of the fuselage surface pressure coefficients is shown to diffuse fairly quickly aft of the sting (see 
figures 4 and 31). Interference of the sting on the predicted total force and moment coefficients results in 
an increase in lift over the free air configuration, an increment that decreases with increasing 
angle-of-attack (figure 33, table 3). The presence of the blade sting also causes an increase in drag at all 
angles-of-attack but 0’ (figure 33) where there is a slight drag decrement. The effect of the sting on 
predicted drag coefficient is shown to increase with increasing angle-of-attack. A stabilizing pitching 
moment increment is predicted as a result of sting interference (figure 33). It should be noted in the charts 
presented in figures 31 and 33 that the data plotted for 8’ angle-of-attack were obtained from solutions 
utilizing flux vector splitting and first order interpolations as opposed to flux difference splitting and higher 
order interpolation. This would tend to introduce more diffusion into the numerical scheme which could 
result in the sharp decrease in lift coefficient as is evident in figure 33. 

A more detailed analysis of the pressure coefficient increment, at 0.95 Mach number and a=4\ 
resulting from the blade sting is presented in figures 34 - 36. For this analysis, only the outboard wing flaps 
were deflected - 10‘ down along the leading-edge and 3' down along the trailing-edge. The increments 
were calculated by interpolating the solution for the free air configuration onto the unstructured mesh 
developed about the configuration with the blade sting modeled. Increments are computed by subtracting 
the interpolated free air levels from the levels of the configuration with the blade sting. For clarity, contour 
lines for AC p =0 are not plotted. The pressure coefficient map presented in figures 34 and 35 shows that 
the sting affects the fuselage lower surface pressures as expected -- a relative compression at the 
leading-edge of the sting, a relative expansion around the sting and a relative compression at the 
trailing-edge of the sting. (The sting has an airfoil cross section (see figure 4).) Notice also the tight 
clustering of contour lines along the leading-edge of the horizontal tail. An unexpected result of sting 
interference is the compression that is seen between the nacelles (figures 35 and 36). (The pentagonal 
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and hexagonal contour patterns are not real but are due to numeric anomalies in the interpolation 
algorithm.) The data presented in figures 31 and 32 suggest that similar surface pressure increment 
trends should be expected for the additional Angles-of-attack analyzed. 

Examination of the previous charts yields some explanation of the behavior of the force and moment 
increment curves of figure 33. The sting interference causes more positive pressures over the lower 
surface of the aft portion of the wing as is shown in figure 35. This would result in increased lift as well as 
more nose-down (negative) pitching moment since the moment reference center is forward of this portion 
of the wing (figure 3). Since the prediction method solves the Euler equations and thus is inviscid, any 
drag predicted is purely a function of pressure. Therefore, an increase in lift due to the sting would lead to 
an expectation that drag would also increase. This is generally the case with the exception of the O' 
angle-of-attack solution. Here, where the sting causes the largest lift increment, a drag decrement is 
predicted (figure 33). The data presented in figures 29 and 30 reveals that the minimum drag in the 
experimental data as well as in the computed drag for the free air configuration occurs near 1 ' 
angle-of-attack. A small increase in lift from the 0’ angle-of-attack value moves the data point to the right 
along the drag polar towards the minimum drag point and therefore to a lower drag level. 

The change in angle-of-attack that results from the presence of the blade sting is shown by the data 
plotted in figure 37. This chart shows the amount that the free air configuration angle-of-attack would need 
to be increased so that the lower surface pressures at the 50% semi-span location would match those 
computed for the configuration with the blade sting modeled. (Any change in lower surface pressure 
coefficient would be an indication of a change in angle-of-attack.) It is evident that the effective 
angle-of-attack increase is indeed less than a degree at O' angle-of-attack. Effective up-wash induced by 
the blade sting is also shown to increase with angle-of-attack. 

The effect of the blade sting on configuration forces and moment when the horizontal tail is deflected 
-5* (leading-edge down) is also presented in figure 33. This figure indicates that deflection of the 
horizontal tail causes a non-linear change in the blade interference from the 8 HT =0' configuration. 

The surface pressure increments for 0.30 Mach number and six degrees angle-of-attack are presented 
in figures 38-41 . In this analysis the inboard and outboard leading-edge flaps were deflected 30' and the 
inboard and outboard trailing-edge flaps were deflected 20’. The increments are again calculated as the 
difference between the sting on results and the sting off results. For clarity, contour lines for AC p =0 are not 
plotted. These charts show that the sting had little effect on the lower wing surface pressures. As a result 
of sting/trailing-edge flap deflection interference, a relative compression, a relative expansion and another 
relative compression, all near the leading-edge of the fuselage-sting junction, is evident in figures 40 and 
41 . Notice the smaller effect of the sting on the aft fuselage at the subsonic conditions as compared to the 
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transonic analysis (compare figure 38 with figure 34). The relative compression between the nacelles 
evident in figure 35 for the transonic conditions is not predicted in the subsonic analysis (figure 39). 

Aft Body Closure Corrections 

The contribution of the afterbody on total force and moment coefficients is shown in figure 42. The 
computational increments were determined by integration of the predicted inviscid pressures of the free air 
configuration over only the aft fuselage and tail surfaces. Corresponding equivalent flat plate skin friction 
drag coefficients are also presented. The experimental increments were obtained by differencing the force 
and moments of the complete configuration and a truncated afterbody configuration. The experimental 
data has been corrected for blade interference utilizing the predicted increments shown in figure 33 and 
tabulated in table 3. Also presented in figure 42 is an increment showing the tolerance level for each force 
and moment increment plot. In general, there is good agreement between the Euler results and the 
corrected experimental data. However, there is some significant differences, the nature of which is not 
fully understood. The differences are partially due to the lack of viscous effects in the computational 
method. The uncertainty in the experimentally data as well as variability in the seal around the sting/ 
fuselage intersection are also possible causes of these differences. 

Conclusions and Recommendations 

From this study, it can be concluded that the USM3D unstructured Euler solver can be used to predict 
the aerodynamic trends for the HSR Reference H configuration. Even though numerical difficulties were 
encountered in the solutions at 8’ angle-of-attack, the code can predict both the levels and trends of the 
experimental wing and fuselage surface pressures across the angle-of-attack range analyzed. Since the 
purpose of this study is to be able to predict sting increments, the inability to predict absolute force and 
moment levels is not too significant since the experimental trends are computed well. It can also be 
concluded from the figures presented that the blade sting causes an increase loading over the aft portion 
of the lower surface of the wing. This results in an increase in lift and drag and a more nose down pitching 
moment. 

For future work, it is recommended that an unstructured grid about a full span configuration be 
developed so that predictions of sting interference at yaw angles may be made. This will necessitate that 
the geometry modeled include less detail than that used in the present study since this semi-span grid 
utilized nearly 200 megawords of computer memory to run. The only other recourse would be for the 
USM3D code to be modified to utilize less memory. It is also recommended that Navier-Stokes solutions 
utilizing the semi-span grid be obtained as soon as the viscous version of USM3D and the corresponding 
grid generation software is available. 
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Figure 2. Close-up of lower surface wing grid. 
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Figure 3. Planform view showing location of experimental wing pressure port rows. 



Figure 4. View of the HSR Reference H aft body showing the empennage and blade sting locations. 
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Figure 5. Comparison of predicted and experimental wing[surface pressure coefficients; M=0.95, a=0 
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Figure 8. Comparison of predicted and experimental wina surface pressure coefficients; M=0.95 a=6 
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Figure 9. Comparison of predicted and experimental wing surface pressure coefficients; M=0.95, a=8 
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Figure 10. Comparison of predicted and experimental wing surface pressure coefficients; M=0 95 

a=4\ 5 ht =-5\ 5 le =0710\ 6t E =073'. 
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Figure 11. Comparison of predicted and experimental wing surface pressure coefficients; M=0.95, 

0=6", 5 H t=-5\ 5 le =0710\ 6 te =073\ 
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Rgure 16. Fuselage pressure coefficient comparison at F.S. 65.306; M=0.95, a=6', 8 H t=0\ 5 le =0710 
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Figure 17. Fuselage pressure coefficient comparison a\ F.S. 65.306; M=0.95, 0t=8\ 5ht=0*, 5 
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Fiqure 18. Fuselage pressure coefficient comparison at F.S. 65.306; M=0.95, a=4’, 5 

8 LE =0710 r 8 T e=073*. 
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Figure 19. Fuselage pressure coefficient comparison at F.S. 65.306; M=0.95. CX=6’ 8ut=-5' 

8 le =071 0 , 5 te =0*/3‘. ht 
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Figure 20. Fuselage pressure coefficient comparison at F.S. 65.306; M=0.30 a=6* 5ut=0 

8 le =30730* 1 5 T e=20720'. 
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Fiaure 21 . Fuselaqe pressure coefficient comparison along a line 45" below the horizontal; M=0.95, 

0C=0', 8 ht =0’, 8 le =0710\ 5 te =073\ 


- 0.5 


s 

• N 

X 



60 65 70 75 80 85 


Fuselage Station, inches 

Fiqure 22. Fuselage pressure coefficient comparison along a line 45' below the horizontal; M=0.95, 

0l=2', 5 H t=0’, 5 le =0710', oj E =073'. 
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Figure 23. Fuselage pressure coefficient comparison along a line 45' below the horizontal- M=0 95 

0=4', 8 ht =0', 8 le =0710', 5t E =073\ 
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Figure 24. Fuselage pressure coefficient comparison along a line 45' below the horizontal- M=0 95 

a=6", 5 ht =o\ 8 le =o7io*, §te=073". 
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Fuselage Station, inches 

Figure 25. Fuselage pressure coefficient comparison along a line 45’ below the horizontal; M=0.95 

0C=8\ 8ht = 0*> 8|_e=0710\ Oxe=0’/3’. 


USM3D Euler, a=4.°, Free Air 

USM3D Euler, a=4.°, With Blade Sting Modeled 

Experiment, a=4.04°, Rn=80x10 6 


Fuselage Station, inches 

Figure 26. Fuselage pressure coefficient comparison along a line 45’ below the horizontal; M=0.95, 

a=4\ 5ht=-5’, 8, £=0710’, 6 te =0’/3’. 






USM3D EuJer, a=6.°, Free Air 

USM3D Euler, a=6.°, With Blade Sting Modeled 

Experiment, a=6.13°, Rn=80x10 6 


Fuselage Station, inches 

Figure 27. Fuselage pressure coefficient comparison along a line 45* below the horizontal; M=0 95 

a=6*, 5 ht =-5\ 8 le =0710*, Ot E =073*. 


USM3D Eufer, 0=6.°, Free Air 

USM3D Euler, a=6.°, With Blade Sting Modeled 

Experiment, ot=6.00°, Rn=10.14x10 6 


Fuselage Station, inches 


Figure 28. Fuselage pressure coefficient comparison along a line 45* below the horizontal; M=0.30 

0C-6°, 5 ht =0°, 8|_e=30730\ 8Ve = 20720“. 










5 le =0710° 5t E =073° 5 ht = 0° 
Mach=0.95 

Note: C D ( USM3D )=C D ( USM3D )+C D sf 
C D sf =0.0073, Rn=80x10 6 


▼ Experiment, Rn=80x10® 

O - USM3D Free Air 

-□ USM3D with Blade Sting 


Figure 29. Comparison of predicted and experimental forces and moment. 












Figure 32. Predicted fuselage pressure coefficients at F.S. 65.306; M=0.95, 
8 H t=0\ 5 le =0710\ 8 te =073’ 


29 



ac d 

0.0010 


-0.0030 L. 
10 o 


Ang!e-of-attack (degrees) 


Angle-of-attack (degrees) 



NOTE: Solution at a=8* obtained 
using Flux Vector Spliting 
and First Order Boundary 
Condition interpolation 


s HT =-s° 


Angle-of-attack (degrees) 


Figure 33. Predicted force and moment increments due to blade sting; M=0.95, 5i p=07l0* < 


















ose-up or a 


uselaae surface pressures increments due to blade sting; M=0.30, ot=6* 
8 HT =0 , 8 LE =30730\ 5 te = 20720’, (AC p =0.05). 
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Skin Friction Drag Estimates 

Rn aCq^ 

10x10® 0.0025 
30x10® 0.0021 
80x10® 0.0018 

Note: Add AC D $f to USM3D AC d only 


•— 

— Experiment, Rn= 10.4 5x10 

B 

experiment, Hn=o 0 . 0 oxi U 

C vna iJiha nl D O *1 

T 

txpenmem, h n= t y.ojxi u 

USM3D, Free Air 


Figure 42. Afterbody / tail surface force and moment increments corrected for predicted blade sting 
interference; M=0.95, 8hj=0\ 8le=0710\ 8te=073* 
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Table 1: HSR Reference H Grids Developed 


5 LE(Ib/ 

8 LE(OB) 

Steob/ 

Stecob) 

Srt 

Sting 

Number of 
nodes 

Number of 
cells 

30730’ 

20720' 

0’ 

No Sting 

223,278 

1,221,259 

30730’ 

20720’ 

0’ 

Blade Sting 

223,380 

1,216,331 

0710' 

073’ 

0’ 

No Sting 

115,306 

839,865 

0710’ 

073' 

O’ 

Blade Sting 

200,007 

1,082,028 

0710’ 

073’ 

0’ 

Circular Aft Body/Sting 

172,244 

929,761 

0710’ 

073’ 

-5’ 

No Sting 

147,623 

797,631 

0710* 

073’ 

-5’ 

Blade Sting 

197,576 

1,068,858 

070’ 

070° 

0’ 

No Sting 

95,687 

505,634 

070’ 

070’ 

0’ 

Blade Sting 

111,899 

591,857 


Table 2: HSR Reference H Configurations Analyzed 


5 LE(Ib/ 

5 LE(OB) 

Steob/ 

^TEfOB) 

5ht 

Sting 

Mach 

Number 

Angle-of- 

attack 

30730’ 

20720’ 

0’ 

No Sting 

0.30 

6’ 

30730’ 

20720* 

0’ 

Blade Sting 

0.30 

6’ 

0710’ 

073’ 

O’ 

No Sting 

0.95 

O', 2’, 4’, 6’, 8’ 

0710’ 

073* 

0’ 

Blade Sting 

0.95 

0\2’,4*,6’,8’ 

0710* 

073’ 

0’ 

Circular Aft Body/Sting 

0.95 

O’, 4° 

0710’ 

073° 

-5’ ; 

No Sting 

0.95 

2°, 4’, 6’ 

0710’ 

073’ 

-5’ 

Blade Sting 

0.95 

| 

070’ 

070’ 

0’ 

No Sting 

2.4 


070’ 

070’ 

0’ 

Blade Sting 

2.4 

4’ 
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Table 3: Blade Sting Interference Effects, M=0.95, 5 L e=0710\ 5t E =073\ 


Angle-of-attack 

mam 

AC|_ 

ACq 

AC m 

O’ 

O' 

0.0428 

-0.0015 

-0.00959 

2" 

O' 

0.0361 

0.00065 

-0.00662 

4' 

O' 

0.0333 

0.0014 

-0.00558 

6' 

O' 

0.0318 

0.0020 

-0.00601 

8* 

O' 

0.0236 

0.0041 

-0.00801 

2' 

-5' 

0.0350 

0.00058 

-0.00577 

4* 

-5' 

0.0318 

0.00166 

-0.00433 

6* 

-5' 

0.0314 

0.00284 

-0.00379 
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